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METHOD  FOR  REMOTELY  MEASURING  FLUCTUATIONS 

IN  THE  OPTICAL  INDEX  OF  REFRACTION  OF  A  MEDIUM 

STATEMENT  OF  GOVERNMENT  INTEREST 

[0001]  The  invention  described  herein  may  be  manufactured  and 
used  by  or  for  the  Government  of  the  United  States  of  America  for 
governmental  purposes  without  the  payment  of  any  royalties  thereon 
or  therefor. 


BACKGROUND  OF  THE  INVENTION 

(1)  Field  of  the  Invention 

[0002]  The  present  invention  relates  generally  to  methods  for 
measuring  index  of  refraction,  and  more  particularly  to  a  method 
for  remotely  measuring  fluctuations  in  a  medium' s  index  of 
refraction . 

(2)  Description  of  the  Prior  Art 

[0003]  Measurement  of  the  index  of  refraction  of  a  light 
propagation  medium  is  traditionally  accomplished  by  collecting 
samples  of  the  medium  and  using  instruments  such  as  an  optical 
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interferometer  or  refractometer  to  measure  the  index  of 
refraction.  Another  approach  involves  measurement  of  quantities 
such  as  temperature,  salinity,  and  pressure,  and  then  using  such 
quantities  in  analytical  models  to  calculate  the  index  of 
refraction.  While  various  methods/instruments  (e.g.,  shadowgraphs 
and  wavefront  sensors)  have  been  used  to  indicate  the  existence  of 
index  of  refraction  fluctuations,  these  methods/instruments  do  not 
have  the  ability  to  remotely  sense  index  of  fraction  fluctuations. 


SUMMARY  OF  THE  INVENTION 

[0004]  Accordingly,  it  is  an  object  of  the  present  invention  to 
provide  a  method  for  remotely  measuring  index  of  refraction 
fluctuations . 

[0005]  In  accordance  with  the  present  invention,  a  method  is 

provided  for  remotely  measuring  fluctuations  in  a  medium's  index 
of  refraction.  From  a  first  location,  a  first  beam  of  optical 
energy  is  focused  at  a  focal  plane  located  at  a  second  location  in 
a  medium  of  interest.  As  a  result,  a  second  beam  of  optical 
energy  is  backscattered  towards  the  first  location.  At  the  first 
location,  a  size  of  the  second  beam  is  determined  where  the  size 
is  indicative  of  strength  of  fluctuations  in  the  medium' s  index  of 
refraction. 
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BRIEF  DESCRIPTION  OF  THE  DRAWINGS 

[0006]  Other  objects,  features  and  advantages  of  the  present 
invention  will  become  apparent  upon  reference  to  the  following 
description  of  the  preferred  embodiments  and  to  the  drawings, 
wherein  corresponding  reference  characters  indicate  corresponding 
parts  throughout  the  several  views  of  the  drawings  and  wherein: 
[0007]  FIG.  1  is  a  top-level  schematic  view  of  a  system  used  to 
remotely  measure  index  of  refraction  fluctuations  in  accordance 
with  the  present  invention; 

[0008]  FIG.  2  is  a  schematic  view  of  a  beam  source  and 
backscattering  measurement  system  in  accordance  with  an  embodiment 
of  the  present  invention; 

[0009]  FIG.  3  is  a  diagrammatic  view  of  a  phase  mask  model  of 
propagation  through  a  medium  with  index  of  refraction 
fluctuations ; ''  and 

[0010]  FIG.  4  is  a  block  diagram  illustrating  the  processing 
steps  used  to  model  propagation  of  electric  field  intensity  at 
the  lens  output  of  the  measurement  system  to  the  system's  focal 
plane  in  accordance  with  the  present  invention. 

DESCRIPTION  OF  THE  INVENTION 

[0011]  The  present  invention  is  a  novel  method  for  remotely 

i 

measuring  index  of  refraction  fluctuations  in  a  light  propagation 
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medium.  The  present  invention  takes  advantage  of  some  basic 
principles  governing  the  interaction  of  a  laser  beam  with  a 
medium.  Briefly,  when  a  laser  beam  is  focused  in  a  medium  such  as 
water  or  air,  index  of  refraction  variations  or  fluctuations 
within  the  medium  can  steer  the  beam,  change  the  divergence  of  the 
beam,  and  break  up  the  beam.  Depending  on  the  strength,  scale, 
and  range  of  the  index  of  refraction  fluctuations,  the  size  of  the 
beam  at  the  intended  focal  point  (or  region)  can  be  significantly 
larger  than  the  diffraction  limited  spot  size.  To  measure  the 
size  of  the  beam  at  the  intended  focal  point,  an  image  of  the 
focused  spot  can  be  made  at  the  source  of  the  laser  beam  by  using 
light  that  is  backscattered  from  suspended • particles  and  molecular 
fluctuations  within  the  medium.  The  image  formed  is  further 
increased  in  size  as  the  backscattered  light  propagates  back 
through  the  index  of  refraction  fluctuations  on  the  return  path. 
The  present  invention  uses  this  change  in  beam  size  in  determining 
index  of  refraction  fluctuations  as  will  be  described  below. 

[0012]  In  order  to  capture  just  a  "snapshot"  of  the  region  of 
interest,  a  variety  of  well-known  optical  techniques  can  be 
employed.  For  example,  light  from  the  focal  point  can  be  isolated 
by  using  a  pulsed  laser  beam  whose  "on"  time  is  short  enough  so 
that  its  extent  in  space  is  less  than  or  equal  to  the  depth  of 
focus  range.  By  time  gating  the  received  signal  to  match  the 
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travel  time  of  the  light  to  and  from  the  focal  point,  the  light 
scattered  within  the  focal  volume  is  further  isolated. 

Measurements  of  the  size  of  the  resulting  image  are  processed  in 
accordance  with  the  present  invention  to  produce  a  measure'  of  the 
index  of  refraction  fluctuations  encountered  during  the  forward 
and  reverse  (i.e.,  backscattered)  propagation  of  the  light  beam. 
[0013]  •  Referring  now  to  the  drawings  and  more  particularly  to 

FIG.  1,  a  schematic  view  of  a  system  used  to  remotely  measure 
index  of  refraction  fluctuations  in  accordance  with  the  present 
invention  is  illustrated.  At  the  heart  of  this  system  is  a  beam 
source  and  backscattering  measurement  system  10  used  to  both  focus 
a  light  beam  12  at  a  focal  plane  100  located  a  distance  (or  focal 
range)  R  from  system  10  and  measure  the  size  of  the  backscattered 
beam  image  that  arrives  at  system  10  along  return  beam  path  14. 
Light  beam  12  and  return  beam  path  14  pass  through  a  medium  200 
(e.g.,  water,  air,  etc.)  that  is  likely  to  experience  index  of 
refraction  fluctuations.  Medium  200  can  be  a  single  medium  or 
multiple -mediums  (e.g.,  air  and  water)  with  a  regular  (e.g.,  flat) 
boundary  between  the  two  mediums .  As  would  be  understood  by  one 
of  ordinary  skill  in  the  art,  system  10  can  be  realized  in  a 
number  of  ways  without  departing  from  the  scope  of  the  present 
invention.  By  way  of  non- limiting  example,  one  possible 
embodiment  of  system  10  will  be  described  later  herein. 
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[0014]  For  simplicity,  the  present  invention  will  be 
explained  for  the  case  of  a  single  focal  range.  However,  it  is 
to  be  understood  that  the  method  described  herein  can  be 
extended  for  the  measurement  and  use  of  data  at  multiple  focal 
ranges  where  such  multiple  measurements  could  be  made 
simultaneously  or  sequentially.  Referring  still  to  FIG.  1,  a 
measurement  of  the  beam  size  at  the  focal  range  R  is  an  indirect 
measurement  of  the  strength  of  index  of  refraction  fluctuations 
in  the  propagation  path  between  the  output  lens  of  system  10  and 
focal  plane  100.  Since  observation  of  the  beam  at  focal  plane 
100  by  system  10  requires  backscattering  of  the  light  by 
molecular  and  particulate  scattering  of  the  light  and 
propagation  back  through  the  fluctuations  along  return  beam  path 
14,  the  size  of  the  beam  image  seen  by  system  10  increases 
further.  If  the  outgoing  beam  nearly  fills  the  output  lens 
aperture  of  system  10,  the  irradiance  seen  by  system  10 's 
detector  can  be  approximated  by  autoconvolving  the  irradiance  in 
focal  plane  100.  In  any  case,  the  size  of  the  detected  beam 
represents  an  indirect  measurement  of  the  strength  of  index  of 
refraction  fluctuations  in  the  path  between  the  output  lens  of 
system  10  and  focal  plane  100.  System  10  provides  the  size  of 
the  detected  beam  to  a  processor  16  where  the  present 
invention's  methods  are  used  to  relate  the  size  of  the  image  to 
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the  strength,  location  in  range,  and  size  of  the  index  of 
refraction  fluctuations. 

[0015]  As  the  strength  of  the  fluctuations  increases,  the 
size  of  the  detected  beam  increases.  If  l(x,y)  is  the  irradiance 
measured  in  the  detector  plane  with  rectangular  coordinates  (r,  y) , 
a  measure  of  beam  size  can  be  the  mean  radius  of  the  beam  given 
by 


r  = 


|  jdxdyl(; 

c,yyj(x-x) 

f+b-yf. 

1 

jdxdyl(x,  v) 

where 


(l) 


_  [  \dxdyl{x,y)x  _  j  \dxdy  I  (x,y)y 

“  -  J  J  and  >’  =  ■ 


x  = 


\\dxdyl{x,y) 


|j dxdyl(x,y ) 


Another  possible  measure  of  beam  size  is  the  mean  square  radius 
of  the  beam  given  by 


r  = 


j"  ^dxdy  I  {. 

x>y%x-* 

M 

iy-yf) 

J 

jdxdy  I  (x,  y) 

(2) 


When  appropriate,  r  and  r'can  be  averaged  over  more  than  one  laser 
pulse . 

[0016]  If  beam  12  is  a  laser  beam  having  a  TEM00  transverse 
mode  structure,  electric  field  intensity  at  the  lens  output  of 
system  10  is  given  by 
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(3) 


-2 

kr 


E_(?,o)^E0e-^  e~J2R 


where 


j  is  the  unit  imaginary  number,  j  =  \— T  ; 

k  is  the  wave  number,  k-  —  \ 

X 

X  is  the  light  wave  length  in  the  medium; 

E0  is  a  complex  constant; 

r  =  xx+  yy  ; 

randy  are  orthogonal  unit  vectors  in  the  (r, y)  plane  of  the  lens; 
R  is  the  distance  to  the  range  gated  focal  plane  100;  and 


Vz 


is  the  e  1  radius  of  the  beam  irradiance  /! 


(ho)-  e_(t,o] 


at 


the  lens  output  of  system  10. 

[0017]  If  medium  200  is  free  of  index  of  refraction 
fluctuations,  the  electric  field  intensity  E_  will  propagate  to 
form  a  diffraction  limited  focused  spot  in  the  focal  plane  at 

the  distance  R  from  the  lens.  In  this  case,  the  irradiance  e 
radius  of  the  Gaussian  shaped  beam  in  the  focal  plane  is  given 
by 


-i 


r_,  = 


Rfi 


a* 


(4) 
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However,  in  the  presence  of  index  of  refraction  fluctuations, 

A n[r,z),  along  the  propagation  direction,  z  ,  the  size  of  the  beam 
will  increase  in  focal  plane  100  or  at  z=R. 

[0018]  Prior  to  describing  the  processing  approach  of  the 
present  invention,,  an  exemplary  embodiment  of  system  10  will  be 
described  with  reference  to  FIG.  2.  A  computer/controller  20 
provides  control  signals  to  a  trigger  22  that,  in  turn,  generates 
an  electronic  pulse  used  to  trigger  operation  of  a  laser  24 .  For 
the  illustrated  embodiment,  it  will  be  assumed  that  laser  24 
generates  a  linearly  polarized  light  pulse  with  a  single  TEM00 
spatial  mode  with  good  beam  quality.  In  accordance  with  an 
additional  control  signal  generated  by  computer  20,  zoom  optics  26 
expand  the  diameter  and  adjust  the  convergence  of  the  laser  light 
beam  to  produce  a  focused  spot  of  the  desired  diffraction  limited 
size  at  a  selected  range  R.  The  diameter  of  the  light  beam  is 
selected  to  achieve  the  desired  index  of  refraction  measurement 
sensitivity  over  the  selected  range.  The  transmitted  light  beam 
further  passes  through  a  polarizing  beamsplitter  28  that  is 
configured  to  allow  passage  of  the  light  beam.  A  quarter 
waveplate  30  converts  the  light  into  circularly  polarized  light. 

A  lens  32  then  focuses  the  transmitted  light  beam  (referenced  by 
arrowheads  12)  through  medium  200  and  onto  particle  and  molecular 
density  scatterers  contained  with  the  focal  region  102  of  beam  12 . 
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Light  that  is  scattered  backward  (referenced  by  arrowheads  14)  is 
then . collected  by  lens  32.  To  the  extent  that  the  polarization  of 
the  light  is  preserved  in  the  scattering  process,  quarter 
waveplate  30  converts  the  returning  backscattered  light  to 
linearly  polarized  light  with  a  polarization  that  is  rotated  90° 
relative  to  the  state  of  the  polarization  generated  by  laser  24 . 
Backscattered  light  with  the  90°  polarization  rotation  is 
deflected  by  beamsplitter  28  into  imaging  optics  34.  Imaging 
optics  34  are  configured  by  computer  20  to  image  the  light  through 
a  shutter  36  and  onto  a  detector  array  38  with  the  desired 
magnification.  A  delay  40  controlled  by  computer  20  turns  on 
shutter  36  to  allow  the  passage  of  light  that  has  traveled  the 
requisite  time  to  and  from  focal  region  102 .  The  open  time  for 
shutter  36  is  set  to  pass  light  that  is  scattered  within  the  focal 
range  of  region  102  and  to  block  light  that  is  scattered  outside 
of  the  focal  range.  The  image  formed  on  detector  array  38  is  then 
provided  to  computer  20  where  range  and  beam  diameter  are 
determined.  Note  that  a  number  of  elements  of  system  10  can  be 
changed  or  modified  without  departing  from  the  scope  of  the 
present  invention.  For  example,  the  use  of  a  broadband  detector 
array  and  time  gating  could  replace  shutter  36  and  detector  array 
38  . 
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[0019] 


Existing  index  of  refraction  fluctuation  studies 


illustrate  the  complex  relationship  between  the  irradiance  I\ r,R) 
and  the  spatial  spectrum  and  strength  of  the  fluctuations.  See, 
for  example,  Andrews  et  al . ,  Laser  Beam  Propagation  through 
Random  Media  (SPIE  PRESS,  1998) .  One  approach  for  determining 

/(r,i?)  in  accordance  with  the  present  invention  uses  a  computer- 
implemented  process  described  below.  The  process  assumes  a 
Kolmogorov  spatial  spectrum  (see  Andrews  et  al . )  for  the  index 
of  refraction  fluctuations.  An  estimate  is  made  of  the  inner  (/0) 
and  outer  (L0)  scales  of  the  fluctuations. 

[0020]  In  terms  of  the  system  embodiment  illustrated  in  FIG. 

2,  the  propagation  path  between  lens  32  and  focal  plane  100  can 
be  broken  into  constant  index  of  refraction  slabs  of  length 
L.  =  Lg  .  As  described  by  Andrews  et  al . ,  a  phase  mask  with 
amplitude  transmittance 


(5) 


where 


is  used  at  the  beginning  of  each  slab  to  account  for  the  index 
of  refraction  fluctuations  within  the  slab.  This  situation  is 
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depicted  in  FIG.  3  for  three  slabs  at  z=0,Lz,2Lz  with  the  input 
field  being  denoted  above  the  arrow  pointing  to  each  slab,  the 
phase  mask  being  denoted  above  each  slab,  and  the  output  field 
being  denoted  above  the  arrow  pointing  away  from  each  slab. 

[0021]  Monte  Carlo  realizations  of  the  0m(r)  phases  can  be 

accomplished  by  using  a  Gaussian  random  number  generator  to 
create  two-dimensional  white  noise.  A  two-dimensional  Fourier 
transform  can  be  used  to  convert  the  white  noise  to  wavenumber 
space  where  the  spatial  spectrum  is  multiplied  by  a  Kolmogorv 
spectrum  before  an  inverse  transform  is  used  to  return  to 

position  space.  An  index  of  refraction  structure  constant  C2  is 
selected  and  the  amplitude  of  the  phase  mask  noise  is  adjusted 
until  the  variance  of  the  phase  mask  fluctuations  satisfies 

8 

=  0.275C2Z|  (6) 

where  CTf  is  the  mean  square  phase.  The  process  described  above 

can  be  used  to  generate  independent  realizations  0m(r)  °f  the 
phase  masks . 

[0022]  With  the  phase  mask  amplitude  transmittance  determined 
as  just  described,  the  process  used  to  model  propagation  of  the 
field  at  the  lens  output  to  the  focal  plane  is  depicted  in  FIG. 

4  for  the  first  slab.  As  shown  in  FIG.  4,  the  electric  field 
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intensity  E_\r, 0]  at  the  output  of  lens  32  is  multiplied  at  40  by 
the  transmittance  ro(r)  of  the  first  phase  mask  42.  The  result 
E+  (r,o)  undergoes  a  two  dimensional  Fourier  transform  at  44  to 
wavenumber  space  to  generate  the  plane  wave  amplitudes  £+(k",o), 

where  K  is  the  projection  of  the  wave  vector  k  =  >c + z  in 
the  (.x, y)  plane  and  z  is  a  unit  vector  perpendicular  to  the  (;t, y) 
plane.  The  transmission  function 


fz(/c)=eiLl'f^ ,  (7) 

corresponding  to  the  plane  wave  phase  shifts  46  multiplies  E+[k, o) 

at  4  8  to  produce  the  wave  number  spectrum  E_[k, Lz)  of  the  electric 
field  at  the  end  of  the  first  slab.  A  two-dimensional  inverse 
transform  at  50  is  used  to  generate  the  electric  field  intensity 

E_[r,L. )  at  the  end  of  the  first  slab.  This  process  is  repeated 

for  each  slab  until  the  electric  field  intensity,  e{t,r)  ,  in  the 
lens  focal  plane  is  determined.  The  relative  irradiance 


/ (r, r)=  |is(rt Rt  is  then  determined. 


[0023]  An  autocovolution  of  /  (r, 7?)  is  performed  to  determine 
the  range  gated  irradiance 


1 DS 


(r,^?)=  Jiiro  l(ro,R.)l(r- ro,R.) 
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(8) 


IDS[r,R)  can  then  be  used  in  Eq.  1  to  determine  the  simulated  mean 

radius  rs  .  The  mean  radius  is  then  compared  with  the 
experimentally  determined  mean  radius.  The  structure  constant, 
C2 ,  is  changed  accordingly  and  the  propagation  process  is 
repeated  until  the  measured  and  simulated  mean  radii  are  in  good 
agreement.  The  final  value  of  C2  is  used  in  Eq.  5  to  determine 

the  mean  square  phase  <J#  which  constitutes  another  measurement 

of  the  index  of  refraction  fluctuations  in  the  propagation  path. 
If  desired,  the  process  can  be  repeated  to  provide  for  ensemble 
averaging  of  the  simulated  beam  size. 

[0024]  The  above-described  process  can  also  be  repeated  for  a 
longer  focal  range.  The  simulation  process  is  repeated  using 
the  structure  constant  determined  to  characterize  the  first 
focal  range  and  as  a  starting  point  for  modeling  the  region 
between  the  two  focal  ranges.  The  search  process  is  repeated 

for  a  range  between  the  focal  planes  until  Cn  and  <7^  are 

determined  for  the  region  between  the  two  focal  planes.  This 
process  can  be  repeated  for  any  number  of  ranges  that  can  be 
achieved,  within  the  limits  imposed  by  attenuation  of  the  light 
within  the  propagation  medium  and  the  laser  power  level,  to  form 
C2  and  <7‘  range  images  of  the  strength  of  the  index  of 
refraction  fluctuations. 
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[0025]  It  will  be  understood  that  many  additional  changes  in 
the  details,  materials,  steps  and  arrangement  of  parts,  which  have 
been  herein  described  and  illustrated  in  order  to  explain  the 
nature  of  the  invention,  may  be  made  by  those  skilled  in  the  art 
within  the  principle  and  scope  of  the  invention  as  expressed  in 
the  appended  claims. 
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IN  THE  OPTICAL  INDEX  OF  REFRACTION  OF  A  MEDIUM 


ABSTRACT  OF  THE  DISCLOSURE 

A  method  is  provided  for  remotely  measuring ' index  of 
refraction  fluctuations.  From  a  first  location,  an  optical  beam 
is  focused  at  a  focal  plane  located  at  a  second  location  in  a 
medium  of  interest.  As  a  result,  a  beam  of  energy  is 
backscattered  towards  the  first  location.  At  the  first  location, 
a  size  of  the  backscattered  beam  is  determined  where  the  size  is 
indicative  of  strength  of  fluctuations  in  the  medium's  index  of 


refraction . 
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FIG.  1 


FIG.  2 
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FIG.  4 


